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Abstract. - Deploying a recently developed semiclassical theory of quasiparticles in the 
superconducting state we study the de Haas-van Alphen effect. We find that the oscillations 
have the same frequency as in the normal state but their amplitude is reduced. We find an 
analytic formulae for this damping which is due to tunnelling between semiclassical quasiparticle 
orbits comprising both particle-like and hole-like segments. The quantitative predictions of the 
theory are consistent with the available data. 



The revival of interest in studying the de Haas-van Alphen (dHvA) effect in the su- 
perconducting state is driven by the hope that this would provide new k- vector dependent 
information about the superconducting gap A(k). Evidently this would be of particular im- 
portance in connection with anisotropic superconductors where A(k) can have lines of zero's 
on the Fermi surface ||. Unfortunately at this stage there is no concensus concerning the 
mechanism of how the experimentally observed oscillations of the diamagnetic response of a 
Type II superconductor come about §-|l|]. Using our recently developed very general semi- 
classical theory of quasiparticles in the superconducting state [jl(| , in what follows we develop 
a semiclassical picture of Landau like orbits of quasiparticles suggested by the simple model 
calculation of Miller and Gyorffy Q. Clearly the long term aim of a semiclassical theory 
is to provide an analogue of the Lifshitz-Kosevich formulae for superconductors. Hopefully 
such a formulae would allow the interpretation of experiments in terms of the Fermi surface 
and the variation of A(k) on the Fermi surface. At this stage we only deal with conven- 
tional superconductors with the usual s-wave pairing. As it happens this is the only case 
for which reliable data already exists. Within the limits of a number of simplifying assump- 
tions the above semiclassical theory provides for Bohr-Sommerfeld like quantisation rules for 
quasiparticles. In particular it allows for the analogue of magnetic breakdown which involves 
tunnelling between distinct semiclassical orbits. We show that there are orbits which enclose 
areas that are precisely the same size as the Landau orbits in the normal state, but involve 
such tunnelling. As will be seen the consequence of these tunnelling events is a damping factor 
in the Lifshitz-Kosevich formulae in agreement with experiments llJ-0] . 

The theory we wish to use seeks the semiclassical spectrum of the following Bogoliubov-de 
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Gennes (BdG) equations 



|A(r)|e-**W (p - eA(r)) 2 - V(r) + e, A ^(r) ^ A V ^(r) 



^(p- f Alrir + nrl-r, A(r)|r— " \ ( „ N (r) ) ( „ A (,-) 

(1) 

where u\(r) and V\(v) are the probability amplitudes that an elementary excitation is a 
quasiparticle and quasihole respectively, A(r) = |A|e l ^ r ) is the order parameter and the 
other symbols have the conventional meaning. For Type II superconductors in large magnetic 
fields A(r) takes the form of the Abrikosov flux lattice |l7[, comprising a periodic array of 
vortices. The effective classical mechanics, Hamilton- Jacobi equations for the quasiparticles, 
corresponding to (|l]) , is described by the following effective Hamiltonians jl6| : 




(P. rl = p • v.lr! + n \: ( £- + \mv*{r) + V{v) - e F ) + |A(r)|», (2) 

where a = ± and mv s (r) = i?iV0(r) + eA(r) is the superfluid velocity. On a constant energy 
surface these determine, implicitly, the functions p Q (r) to be used in the Bohr-Sommerfeld 
quantisation rule. To simplify matters we take the crystal potential V(r) to be a constant. It 
is clear from the single vortex solution |l6j, that the detailed shape of |A(r)| is not essential. 
This fact will be used below. By contrast the role of the line of phase singularities [|l6| , which 
runs along the vortex core, needs more careful consideration. The topologically non-trivial 
behaviour of the phase gradient • dr ^ for any path enclosing lines of phase singu- 

larities) causes the phase gradient, V0, to cancel out (on the average) the increase in A(r) 
across the cell. Consequently the Abrikosov form for v s (r) is periodic. The Hamiltonian, (g), 
then describes quasiparticles which correspond to spinor Bloch waves . On the other hand 
the topologically non-trivial behaviour of A(r) arises from the superposition of topologically 
trivial solutions to the Ginzburg-Landau theory, so long as the coefficients are suitably deter- 
mined from minimization of the free energy for the nonlinear problem |]l7f . We can therefore 
consider the solution of the BdG equations for a topologically trivial A(r) = |A(y)| exp(ik x x), 
with \A(y)\ = A exp(-i(y - fff) 2 /( 27r )~ lrf2 ) and k x = 2-nnjd (d the flux lattice cell size, 
n an integer). A full Bloch wave solution to equation (]l|) can then be constructed by form- 
ing an appropriate superposition and determining the coefficients from self-consistency. The 
single- valuedness of each solution ensures the resulting superposition is also single- valued, and 
the former requires the quantisation of a Bohr-Sommerfeld integral for quasiparticles with the 
A(r) given above. 

Taking this approach, and working in the Landau gauge, we rewrite (0) as 




£"'lp.r) = /V_i/;H-n\; ( £L- + l mVsx (y)-e F S j +\A(y)\*, (3) 

where mv StX (y) — hk x /2 + eA x (y) and we have approximated B(r) = V x A(r) by a constant 
(averaged) B-field. The quasiparticle dynamics described by (||) are integrable and we can 
deploy the full semiclassical machinery developed in p6| . This we now do. 

The momentum branches defined by the Hamilton-Jacobi equations E a (p, r) = E a are 

Pp a (y) = \jpl- P 2 X -Vl- m^ c (ft - yf ± 2m^J (E° + P x co c (y - - ^(^ 



(4) 
nd 

(n integer). Less obviously, P x , also behaves like an orbit centre (see below). It is natural 



We see immediately that |^g- behaves like an orbit centre, though restricted to = nd 
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Fig. 1 - Phase space orbits defined by p^' a given in (Q) with |A(y)| = |A|. The solid line is for a 
particle-like excitation, the dashed line is for a hole-like one. Arrows indicate the direction of the 
velocity, (a) P x = 0. The turning points along the y-axis are given by p^' a = 0. (b) P x 7^ and 
small. The orbits are shifted in opposite directions. The vertical dotted lines indicate the position of 
the Andreev- like turning points, y+(P x ) &n&y-(P x ), see equation (^|). (c) As P x increases further the 
Andreev-like turning point y+, for which p^' a 7^ 0, is reached before the normal turning points given 
by p^' a — 0. MAS takes place (see text), (d) P x increases even further. Now MAS also occurs at y_ 
creating a second particle-hole orbit. The two orbits are separated in real space by Sy — 2jA|/ P x u) c . 



to restrict it to one flux cell, i.e., — | < < +4. Then, since lies within a cell, and 

' '2 — eB — 2 ' eB 5 

represents translations by lattice cell vectors, the pair combined allow us to centre orbits 
anywhere in the sample. As we will see, different placements of orbits within a cell lead to 
different E\, i.e., degeneracy within a cell is lifted. On the other hand translations of such 
an orbit by lattice vectors (different choices of k x ) leave E\ = E\(P X ) invariant. One further 
simplification replaces |A(y)| by (|A(y)|) = |A| in the region of interest. Outside this region 
we expect to set |A(j/)| = but, since from (0), any choice in the range < A(y) < |A| for 
y — Tik x /2eB > d makes little difference, we can for simplicity still take |A(y)| = |A|. Taking 
all these observations into account the constant energy orbits corresponding to these relations 
are the same as those studied by Burmistrov and Dubovskii [^3|. Firstly we investigate the 
orbits and spectrum for P x = 0. Fig. [l](a) shows typical classical phase space orbits defined 
by (Q) with P x = 0. The corresponding Bohr-Sommerfeld quantisation condition yields 

E±(p,) = ±^ (hu c (n + + || - CJ p) 2 + | A|* . (5) 

This is the Landau level spectrum, shifted into two square root singularities (Fig. ^), which 
was studied previously by Miller and Gyorffy || and Miyake [§). As it turns out this theory 
is over simplified but nevertheless it contains the basic physics. The picture is one of Landau 
levels which march across the gap || giving the same frequency as in the normal state, but 
with an extra damping, due to |A| ||||, given by R sc = aKi(a), with a = 2n£j^-, (£ integer), 
where Ki is the Bessel function for imaginary argument. Whilst R sc (a) fits the experimental 
data it predicts too much damping for realistic values of |A|. However, as was pointed out 
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in ref. |E| P x ^ changes the spectra drastically because, as is clear from (g) 'Landau level' 
energies depend upon P x . The rest of this paper deals with these complications. 

Thus we consider the orbits and spectrum for P x ^= 0. Figure 0(b), (c) and (d) demon- 
strate the marked change in the structure of the phase space orbits as P x increases from zero 
(figure 0(a)). To understand the behaviour of the branches of P^' a , it is helpful to view the 
orbit structure as resulting from the competition between two different types of turning point. 
The usual turning points, P^' a = 0, are due to the harmonic oscillator confining potential, 
^muj^y 2 , (y — y — hk x /2eB) in (0), and are those seen in figure (a) and (b). However the 
presence of the second square root in (0), y/{E + P x uj c y) 2 — | A| 2 , which we denote by e{y) 
provides an alternative reflection mechanism for which p^' a ^ 0. If e(y) becomes zero before 
the normal turning points are reached the momentum becomes complex and the wave function 
will correspondingly exhibit evanescent decay. Classically the excitation undergoes reflection. 
These new turning points given by e(y) = are 

±\A\-E 

This reflection process is analogous to Andreev scat tering [^9[. Andreev scattering is due 
to scattering from inhomogeneities in |A(r)| Thus \JE 2 — | A(r)| 2 — -> as |A(r)| — > E from 
below. The reflection process we consider has |A| = constant, but due to the non-zero vector 
potential the extra term, P x ui c y, can still cause e(y) — > 0. (It is for this reason that the 
detailed shape of |A(y)| is not essential.) To emphasise the similarity of this process with 
Andreev reflection, and the role of the vector potential we shall call the scattering mechanism 
Magnetic Andreev Scattering (MAS), and y± Andreev-like turning points. 

The appearance of MAS when P x is varied (Fig. 0(c) and (d)) can then be explained as 
follows. From (0) we see that as P x — > 0, y± — > — oo (E > |A|). When P x becomes non-zero 
the Andreev-like turning points move in from — oo towards the normal turning points. We also 
see, as stated above, that P x acts essentially like an orbit centre, the particle-like and hole-like 
orbits (see below) being pulled in opposite directions due to their opposite charge (Fig. 0(b)). 
The orbits first intersect when Py{y+) = Py{y+) i-e., when the Andreev-like turning point 
reaches the normal turning point. For a more physical picture of the nature of these turning 
points remember that the spinor (the eigenfunction to (0)) being transported along a given 
trajectory has both particle, and hole, |?;|, amplitudes and thus carries an effective charge 
e* = e(|w(r)| 2 — |w(r)| 2 ). Evidently e*(r) is a function of position and in particular changes 
from e* > along p+ (Fig. 0(c)), through zero (at the turning point) to e* < along p~ . So 
the effective charge of the excitation changes sign at ?/+ and correspondingly the direction of 
circulation in the magnetic field is reversed. A particle-like excitation is reflected as a hole-like 
excitation, and a given orbit has both particle-like and hole-like segments. 

Having discussed the appearance and interpretation of the MAS orbits we now turn to 
their quantisation. The spectrum obtained from semiclassical quantisation is shown in Fig. 0. 
Firstly note that at P x = there are Landau levels corresponding to (0) which are excluded 
from the gap. As we increase P x we reach the situation shown in Fig. 0(c). The sensitive 
dependence of the turning point y+(P x ) upon P x results in a spectrum which is highly P x 
dependent. Increasing P x further results in the situation shown in Fig. 0(d). The dramatic 
new feature of the spectrum for P x ^ is the existence of states inside the gap. These states 
are new features of the theory which were not there in the P x = case. However despite the 
existence of these states it is clear that they cannot on their own account for the experimental 
facts. These states have considerably different phase space (and hence momentum space) 
area, so that the dHvA frequency for each individual orbit, which is related to this area, will 
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Fig. 2 - Landau level spectrum denned by (g). The BCS density of states is included to emphasise 
how it is broken upin a magnetic field into Landau levels, which are nonetheless pushed apart by the 
gap, |A|. 

Fig. 3 - Semiclassical spectrum obtained by quantising the orbits in Fig. [I]. Region A corresponds to 
quantising the orbits in Fig. |l](a) and (b), whilst region B and C correspond to quantising the one 
or two scattering orbits in Fig. |l|(c) and (d). The spectrum first becomes gap less when P x = p m in ~ 
\A(B)\/v F . 



be very different from that in the normal state. To explain the experiments we in fact require 
one further bit of physics. That is a phenomenon analogous to 'Magnetic Breakdown', pl|,p2| 
which we shall now discuss. 

We observed that the separation of the two orbits in Fig. 0(d) is given by Sy — y+ — y_ = 
2|A| / P x cu c and it thus inversely proportional to both P x and the B-field (through oj c = eB/m). 
For a fixed B, as P x increases this separation in real space rapidly shrinks and quasiparticle 
tunnelling becomes a viable option. Purely on physical grounds we sec that tunnelling be- 
tween orbit segments (Fig. ^) can reproduce an orbit in phase space with approximately the 
normal area. To quantise explicitly the tunnelling orbit the Bohr-Sommerfeld integral can be 
represented as a complex contour integral. The Riemann surface for the problem is somewhat 
similar to that of Fortin et al [^(| . In our case the MAS orbits live on two different sheets of 
the Riemann surface but the tunnelling orbit remains on only one sheet. Consequently the 
integral is equal to the residue at oo, on this sheet, and we find 

E n (p z ) = hu c (n+^J +J^-£f- (7) 

The spectrum is P^-independent and the tunnelling orbits have exactly the normal state 
frequency. However the tunnelling coefficient, which depends upon the imaginary part (Im) 
of the momentum, is 

T(P x )=e Jy -^ y , (8) 

and clearly depends upon P x and B. For fixed P x , Sy cx 1/B — ► as B increases, so that 
the possibility of tunnelling is enhanced. We can therefore view the tunnelling as 'Magnetic 
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Fig. 4 Fig. 5 



Fig. 4 - The tunnelling orbit (solid line) comprised of segments from both of the MAS orbits. Tun- 
nelling from y+ to y- results in an orbit with precisely the same phase space area as for an orbit in 
the normal state. The spectrum quantising such an orbit is given in equation ([?]). 

Fig. 5 - Fit of the current theory for the extra damping, R sc , in the superconducting state to the 
data for VaSi using |A| = 4.4 meV. The upper critical field is B C 2 = 18.5T. 



Breakdown of MAS orbits'. Magnetic Breakdown in the normal state is a well studied phe- 
nomena [^0|. The oscillatory magnetisation formulae is modified due to the tunnelling. To 
account for this we must include the factor T 2k (every revolution of the orbit we pick up 
2 tunnelling coefficients, and in general we have k revolutions). The P x dependence of (||) 
requires the J2p T 2k to be carried out when calculating the magnetisation. (An additional 
J2k yields the usual degeneracy factor in the formulae for the magnetisation.) Observing 

that Px.max ^ Pf we nn d / dy'lm Py{y') ~ f ^r^p"; where vp is the Fermi velocity. 
Using |A(£?)| = A(0)(1 — B/B^) 1 ^ 2 — > as B — > B C 2 we have the important observation that 
T{P X ) — > 1 for all P x as B — > B C 2. Thus the additional damping vanishes as |A| — > as is to 
be expected. Note, tunnelling is not a small correction, it is zeroth order. The normal state 
is recovered for T(P X ) — > 1. 

The main contribution to ^ p T 2k (P x ) is obtained for the largest value of P x = eBd/2 
and is given by 

T 2k (P x = eBd/2) = e - 2i ^ fc ££^, (9) 

where A = (2eBh)^ 1 ^ 2 . Interestingly this result is very similar to that of Wasserman and 
Springford More generally, the extra damping factor in the superconducting state is 

found to be 

1 fPd(B) _2k 

R SC (B) = - / dP x e h 

Pd Jpmin{B) 

where pd — eBd/2, and p m in(B) = \A(B)\/vf is the smallest value of P x for which the 
spectrum is gapless (see Fig. |||). Although the above results have been derived for a free 
electron model with short range attraction, constant (averaged) pairing potential (|A|) and 
uniform magnetic field B, provisionally the formulae for the damping factor given in ( fL0| ) can 
be compared with that deduced from experiments. We replace vf by its orbitally averaged 
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quantity, and use (1/vp) — ra^j \j2heF with m& = 0.9m e , F = 1560T, and |A| = 4.4 meV to 
compare ( |To| ) with the data for VaSi in figure [|. Even at this early stage we have a satisfactory 
fit to the data. Clearly, |A| in the above analysis is an effective quantity whose value cannot 
be expected to agree with the zero field |A| of 2.6 meV (see for example ref. 0). A more 
detailed analysis of the existing data as well as a more complete presentation of the theory 
will be given elsewhere F). 

To summarise we developed a semiclassical theory in the superconducting state. We con- 
sider new semiclassical orbits involving 'magnetic breakdown' of MAS orbits. These new 
orbits, when quantised, have precisely the normal state frequency. The tunnelling involved 
gives rise to an extra damping in the superconducting state in agreement with experiments in 
such classic superconductors as V3SL When generalised to include period crystal potential, 
a more complete description of the flux lattice, and anisotropic pairing Q , the above theory 
will also be suitable for analysing experimental results on exotic superconductors when ever 
these become available. 

* * * 
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